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Abstract 



An effective field theory for the three-body system with large scattering 
length a is applied to three-body recombination into deep bound states in 
a Bose gas. The recombination constant a is calculated to first order in 
the short-distance interactions that allow the recombination. For a < 0, 
the dimensionless combination ma/ '(ha 4 ) is a periodic function of In |a| that 
exhibits resonances at values of a that differ by multiplicative factors of 22.7. 
This dramatic behavior should be observable near a Feshbach resonance when 
a becomes large and negative. 



One important factor limiting the achievable density in Bose-Einstein condensates 
(BEC's) of trapped atoms is the loss of atoms through 3-body recombination [|l],0. Such 
losses occur when three atoms scatter to form a molecular bound state and a third atom. 
The kinetic energy of the final-state particles allows them to escape from the trapping po- 
tential. This process provides a unique window on 3-body dynamics of cold atoms. In a 
BEC, 3-body loss rates could also reveal new collective mechanisms involving a molecular 
condensate ||. 

The recombination event rate can be parametrized as u rec = an 3 , where a is the recombi- 
nation constant and n the density of the gas. The calculation of a in general is a complicated 
problem, because it is sensitive to the detailed behavior of the interaction potential. If the 
s-wave scattering length a is large compared to the size of the atoms, however, the problem 
simplifies. Assuming that a is the only important length scale, dimensional analysis implies 
a = Cha A /m, where m is the mass of the atoms and C is dimensionless. For a > 0, there 
can be 3-body recombination into a shallow s-wave state with binding energy h 2 / (ma 2 ). For 
either sign of a, there can also be recombination into deep molecular bound states. The 
coefficient a can then be split into two contributions as and ao coming from recombination 
into the shallow and deep bound states, respectively. 

The 3-body recombination into the shallow bound state has been studied in Refs. [fj [7J . 
For the s-wave contribution, C is an oscillatory function of In a: 



where s ~ 1.0064 and A* is the 3-body parameter introduced in Ref. 0. In the adiabatic 
hyperspherical approximation, the oscillatory behavior can be explained by interference 
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between paths connecting two hyperspherical potentials ||||. Within the framework of 
effective field theory, this unique behavior is due to scaling violations created by a 3-body 
interaction required for consistent renormalization J7|. For a < 0, there is no shallow bound 
state (there is a low-energy virtual state instead) and as = 0. 

Effective field theory (EFT) is a powerful method for describing systems composed of 
particles with wave number k much smaller than the inverse of the characteristic range R 
of their interaction. EFT focuses on the universal long-distance aspects of the problem, by 
modelling the interactions as pointlike || . For wave numbers fc < 1/i?, one can expand in 
powers of the small variable kR. Generically, the scattering length a is comparable to R, and 
the expansion is effectively in powers of ka. In the case of large scattering length a 3> R, 
the dependence on ka must be treated nonperturbatively. In the EFT, a single 3-body 
parameter is necessary and sufficient to renormalize the ultraviolet divergences at leading 
order in kR @. The scattering length and this 3-body parameter describe all low-energy 
3-body observables up to errors of order R/a. 

In this Letter, we use EFT methods to study the recombination constant old for deep 
bound states. The deep bound states are outside the range of the effective theory and 
cannot be treated explicitly. Nevertheless, their contribution to 3-body recombination can 
be calculated within the EFT by using the optical theorem (cf. Fig. 0(a)). The 3-body 
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FIG. 1. (a) Illustration of the optical theorem. Two-body bound states are indicated by the 
shaded double line, (b) Diagrams for the contribution of the deep bound states. The single (double) 
lines represent the exact propagators for the field ip (d) and the triangle denotes an insertion of 
ihf. 

recombination channels generate imaginary parts in the 3 — > 3 scattering amplitude. The 
imaginary part from recombination into the shallow bound state is dynamically generated 
in the EFT. The imaginary part from 3-body recombination into deep bound states can be 
taken into account by adding a local interaction term to the EFT. 

A similar approach has been used successfully to describe the effects of the annihilation 
decays of positronium in QED and heavy quarkonium in QCD within the framework of non- 
relativistic EFT's fTTl. The annihilation of positronium into photons involves intermediate 



states with relativistic electrons and positrons that are outside the range of validity of the 
EFT. However, the effects of the annihilation process on e + e~-scattering at low energies can 
be taken into account systematically by adding local 4-fermion operators with imaginary 
coefficients to the effective Lagrangian. In the case of recombination into a molecule with 
binding energy B ^> h 2 /(ma 2 ), the molecule and the recoiling atom will emerge with large 
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momenta of order \JmB which is outside the domain of the EFT. Up to corrections of or- 
der p 2 /(mB), where p is the scale of the momentum of the three atoms, the effect of the 
recombination process on the 3^3 scattering amplitude can be taken into account with 
a local 3-body interaction term whose coefficient has an imaginary part. The interaction is 
local, because the three atoms have to approach to within a distance of order h/ymB <C a 
to recombine. The imaginary 3-body coefficient accounts collectively for recombination into 
all deep bound states. It could in principle be calculated if the 2- and 3-body potentials de- 
scribing the interactions between atoms were known with sufficient accuracy. In practice, it 
has to be determined from experimental data. Once this parameter is fixed, the dependence 
of an on a and A* can be predicted. 

In Ref. [J/J, we have calculated a$ for a large, positive scattering length a. In the following, 
we extend this work to calculate atp for large a of either sign. For simplicity, we now set 
H — 1. We start with a local Lagrangian for a nonrelativistic boson field ip and an auxiliary 
field d with the quantum numbers of two bosons [0,0: 

- -7=(dV^ + + {h + iti) rftdVty + • • • > 

where h! is the free parameter accounting for recombination into the deep states. The dots 
denote terms with more derivatives and/or fields; those with more fields do not contribute 
to the 3-body amplitude, while those with more derivatives are suppressed at low energy. 
The couplings g and h in Eq. @) can be eliminated in favor of a and A*, respectively. The 
atom propagator has the usual nonrelativistic form i/(u — p 2 /2m). The 2-atom scattering 
amplitude for incoming atoms with momenta ±k has the form (— 1/a — ik) -1 . The exact 
propagator for d has a pole at uj — — 1/ (ma 2 ) +p 2 / (4m) corresponding to the shallow bound 
state (7HJ. 

Before describing the calculation of ao, we summarize the calculation of as in Ref. 0. 
The coefficient for recombination into the shallow bound state can be expressed as as = 
5127T 2 \t(pf)\ 2 /(^/3m), where t(p) is the amplitude for the transition between three atoms 
at rest and a final state consisting of an atom and a shallow s-wave state with momentum 
p in their center-of-momentum frame. This amplitude is evaluated on shell at the value 
Pf = 2/(^/3a) prescribed by energy conservation. However, t(p) is also defined at off-shell 
values of p as the solution of an integral equation that contains the 2- and 3-body interaction 
terms 0. Note that t is related to the amplitude T of Ref. via t(p) = mT(p)/(487r 3 / 2 v /a). 
The integral equation is regularized by a momentum cutoff A. All observables can be made 
independent of the cutoff by adjusting the 3-body coefficient h(A) as a function of A. The 
evolution of h with A approaches a limit cycle: A 2 h varies periodically between +oo and 
— oo as A increases by multiplicative factors of exp(7r/so) ~ 22.7 0. Since all observables 
are cutoff independent, we can choose a A for which h(A) = 0. The 3-body parameter A* 



then appears in the upper limit of the integral and we obtain a renormalized equation [11 
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with A n = A*exp[(mr + arctan(l/s ))/<So] an d s o ~ 1-0064. The integer n should be large 
enough that errors of 0(l/(A n a)) can be neglected. In practice, one can often choose n = 1. 
Solving the integral equation for t(p) and evaluating the solution as described above, we 
obtain the result in Eq. ([I]). The integral equation (|3]) also has a solution for a < 0, but 
there is no shallow bound state in this case and as = 0. 

In order to describe the recombination into the deep bound states, we use the optical 
theorem as discussed above and include an imaginary part h! in the coefficient of the 3- 
body contact interaction. We assume that the h! term can be treated as a perturbation 
at first order. However, the external legs have to be dressed with initial and final state 
interactions to all orders. We have to calculate the diagrams shown in Fig. |l|(b), where the 
blob represents the solution of Eq. (||) and the triangle represents an insertion of ih! . This 
leads to 
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where we have defined h'/(mg 2 ) = h\j A 2 n . Eq. @ gives a renormalized expression for an 
up to corrections of 0(l/(A n a)). Since t(q) is known numerically from solving Eq. (|3|), the 
calculation of an entails the numerical evaluation of the integral in Eq. (^). 

There are now two cases to consider: a > and a < 0. We first consider the case a > 0. 
Evaluating Eq. ([|), we find that the amplitude is dominated by a linear divergence in the 
integral and is proportional to A n a. The factors of A n cancel and the result is proportional 
to a 4 with a coefficient that, within the numerical accuracy, is independent of A*a: 



a D = 0.041/ii— (a>0). 
m 
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In Fig. 0(a), we show both as and an for a > 0, with hi adjusted such that an = 10. The 
recombination into the deep bound states fills in the zeroes of as, so that they become at 
best local minima of the recombination rate. 

We now turn to the more interesting case a < 0. Evaluating Eq. (^), we find that the 
amplitude is again dominated by a linear divergence in the integral. However, the coeffi- 
cient of a 4 /m is not a pure number but exhibits scaling violations, an has the remarkable 
dependence on A*|a| shown in Fig. ||(b). The coefficient diverges at a sequence of values 
of A*|a| that are equally spaced on a logarithmic scale. To a high numerical accuracy, the 
results can be described by the formula 
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The oscillatory dependence on In \a\ is in accord with a general scaling law by Efimov 
There are divergences whenever sq ln(A^ |a|) + 1.77 is an integer multiple of ti. They occur 
because a 3-body Efimov state [O has been tuned to threshold by the variation of a. The 



divergences are artifacts of the first order perturbation theory in h\, which treats the Efimov 
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FIG. 2. The recombination coefficients an (solid line) and as (dashed line) as functions of 
|a|A* for a > and a < 0. The free parameter h\ in ar> was chosen such that old = 10 for a > 0. 



states below the 3-body threshold as having sharp energies. However, these Efimov states 
aquire a width that depends on hi from recombination into the deep states This 



will change the divergences into resonances at which the recombination rate is enhanced 
by a factor {Ti 2 /Tma 2 ) 2 , where T is the width of an Efimov state at the 3-body threshold. 
First-order perturbation theory in h\ is justified if T <gC ^ /(ma 2 ), except within T of a 
resonance. The determination of the height of the resonance peaks requires a calculation 
that is nonperturbative in h\. Note that the elastic 3-body scattering amplitude will also 
exhibit these resonances if a < 0. For a > 0, an Efimov state at the 3-body threshold 
has a width of order h 2 /(ma 2 ) from recombination into the shallow 2-body bound state. 
Consequently, there should be no resonances in the recombination rate. 

The EFT results for a depend on two parameters A* and hi that are determined by 
interatomic interactions at distances much smaller than |o|. They have predictive power 
only if a can be varied as a function of some external parameter and if the dependence of 
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A* and h\ on that parameter is known. One case in which there is predictive power is near 
a Feshbach resonance, where a is varied by tuning a molecular bound state close to the 
threshold of two atoms by using an external magnetic field B. Short-distance parameters 
must be smooth functions of B. Over the narrow interval of a Feshbach resonance, they 
can be approximated by a linear or even constant function of B. Linear dependence of the 
short- distance parameter g in fl2|) corresponds to a(B) = an + A/(£> — B ), where an is the 
off-resonant scattering length and A is the width of the Feshbach resonance. The short 
distance parameters A* and hi can be approximated by constants near the resonance at B . 
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FIG. 3. The atom loss rate K 3 = 3a for a gas of ultracold 85 Rb-atoms. The dots (triangles) 
represent data (upper bounds) for K 3 from Ref. 0. The EFT predictions for K 3 s (solid lines) and 
Kzd (dashed lines) are shown for A* = 5/a.u. (black curves) and A* = 20/a.u. (grey curves). The 
arbitrary normalization factor in K^d is set to hi = 24.7. 



In Fig. |H we compare our theory with measurements of the atom loss rate coefficient 
K 3 = 3a in an ultracold gas of 85 Rb-atoms by Roberts et al. ||. In 85 Rb, there is a 
Feshbach resonance with width A = 11.5 G at a magnetic field B = 155.2 G 0. In Fig. |3|, 
the scattering length is negative below the Feshbach resonance at B Q = 155.2 G and positive 
above. Consequently, the shallow bound state contributes to K 3 only above the resonance. 
The free parameter h\ = 24.7 has been fixed to reproduce the overall magnitude of K 3 below 
the resonance. For a > 0, this value of hi leads to a contribution of the deep bound states 
to the loss rate (K^d) that is much smaller than the contribution of the shallow bound state 
(K$g) except near the minima of -^35. One might be able to determine h± experimentally 
from the depth of these minima. In Fig. ^|, we have plotted K 3 £> and -^35 for two values 
of the 3-body parameter: A* = 5/a.u. and A* = 20/a.u.. The parameter A* would be 
determined very accurately if a local minimum or maximum in the loss rate were observed. 

Divergences in the 3-body recombination rate for a < have been observed previously 
by Esry et al. || . They calculated a numerically for various 2-body potentials, most of them 
with a single bound state (either shallow or deep). They varied a from large negative values 
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to large positive values by tuning the depth of the potential. For a > 0, a/a A increases to a 
maximum at some value a max and then decreases to an interference minimum as predicted 
by Eq. (JTJ) . For a < 0, a/a 4 diverges at some value a div . The ratio |cidiv|/ a max ranged 
from 1.6 to 2 for the potentials studied. Since the maxima of as/ a 4 and the divergences of 
afl/a 4 in Fig. ^] appear at the same values of \a\, the EFT predicts |adiv|/a m ax = 1- The 
discrepancy between the results of Ref. and the EFT prediction can be explained by a 
variation in the effective short-distance parameters A* and h\ as the depth of the potential 
is tuned to vary a. Near a Feshbach resonance, the short distance behavior of the potential 
is essentially unchanged as B is varied. Thus A* and h\ should be nearly constant and the 
EFT prediction should hold as long as \a\ is much larger than its off-resonance value. 

We have shown that the EFT approach to the 3-body system with large scattering 
length can be applied to recombination into deep bound states. For a > 0, this contribu- 
tion fills in the interference minima of the recombination rate into shallow bound states, 
transforming them into local minima of the total recombination rate. For a < 0, the EFT 
approach makes the dramatic prediction that the recombination rate for zero-momentum 
atoms should exhibit resonances due to Efimov states at a sequence of values for a that 
differ by multiplicative factors of exp(7r/so) ~ 22.7. Provided the widths of these Efimov 
states are numerically small compared to h 2 /(ma 2 ), the resonances should be observable as 
sharp peaks in the recombination rate near a Feshbach resonance. 
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